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ABSTRACT 


What makes a hump effective? - this is a difficult question. 
In this analysis we compare humps by modelling passenger comfort and 
safety along with the rattle space or wheel well clearance. Design 
considerations for humps involve factors like type of vehicle, 
allowable speed limit, acceptable discomfort level of passenger etc. 
A mathematical model representing the vehicle and driver system is 
developed. The dynamic response of this system is evaluated through 
computer simulation both in time and frequency domains for a class of 
humps. It was found that the dynamic effect of humps vary 
considerably with the type of the vehicle. Size of the hump also 
plays a major role on its dynamic effects. Simulation results show 
that humps cannot be designed that will have little or no effect at 
slow speeds, and the difference in the effects at low and high speeds 
is not very much. This indicates that humps are poor means of 
controlling speed, and that other mechanisms should be sought for 


this purpose. 
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CHAPTER I 


INTRODUCTION 

Road humps are often used for controlling the speed of 
vehicles at sensible spots such as school zones, turning curves, 
residential areas etc. In order to keep vehicle speeds below the 
limiting speed, humps should be designed efficiently. Not much work 
has been carried out in this area to specify guidelines to the road 
designer for designing humps to control traffic. A good hump design 
should have minimal negative effect on the system at lower speeds and 
significant effect at higher speeds. How well can humps meet this 
goal? 


Given a vehicle type, a desired speed limit, and a certain 
discomfort profile, what type of hump should a road designer use? 
This is the question we set out to answer. In order to model humps 
efficiently, one must know the class of vehicles for which the speed 
limits are to be enforced, the level of discomfort that can be 
withstood etc. Typical hump parameters include maximum height, 
maximum span, hump profile and the peak curvature rate. As it is not 
possible to construct and test all types of humps, simulation is done 
to identify the dynamic behavior of the system for travel across 


different types of humps. 
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This work has been divided into four chapters. In Chapter 1 
various hump parameters that are to be considered are discussed. In 
Chapter 2, we discuss the different vehicle models. Our model which 
allows more degrees-of-f reedom ( 8-DOF including driver) than the 
previous works [6] is also discussed. The solution methodology is 
also presented with 1-DOF model as an example. These measures are 
some what subjective e.g. [3, 7 and 14] have different views on 
these. Comfort charts have been derived for sustained vibrations; 
but the non-harmonic shock input of a hump does not fit this 
paradigm. On the other hand, breaking contact occurs quite frequently 
at microscopic scales and thresholds for sensible loss of contact are 
difficult to determine. The criteria for evaluating vehicle 
performance and passenger comfort and simulation results for various 
hump cases are dealt in Chapter 3. Concluding remarks along with 
scope for future work are embodied in Chapter 4. In Appendix A the 
elements of the model equation matrices are given. Two sample vehicle 
parameters representing a light vehicle like a passenger car (Make-1) 
and another medium size vehicle like a van (Make-2), are listed in 
Appendix B. 

1.1 HUMP PARAMETERS 

Dimensions and configurations of hump have a major influence 
on vehicle and driver. Generally a hump is specified by the following 


parameters . 
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1 . 1.1 HUMP PROFILE 


The hump profile refers to the geometry of the hump. Such a 
profile can be circular (i.e., a segment of a circle), parabolic, 
cycloid, harmonic or connected segments of straight lines. Ideally a 
hump profile should at least have smooth entry and exit 
(characterized by zero end conditions i.e., zero height at ends and 
zero slopes at both ends and at mid point), and must be symmetric 
about the center line of its cross-section in order to produce the 
same effects while crossing it from either of its directions. For a 
given profile, the hump parameters are maximum height and span. 

1 . 1.2 OTHER PROFILES 

The following fourth order equation represents half of the 
hump section. It can be adapted according to the requirements of zero 
end conditions and to obtain desirable curvature .The main advantage 
of such a profile is that within the chosen length and height, a 
large variety of shapes can be generated by changing the curvature at 
the ends. The equation of a fourth order spline can be expressed as 

Yi(X)= a.X 4 + b.X 3 + c.X Z + d.X + e 


CD 
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with zero-end conditions i.e., Yi(0) = Yi(0) = Yi(L) = 0, and with 

ft 

Yi(L) = H, Yi( 0 ) = », we can relate the coefficients a, b, c, d and e 
of the above equation with hump parameters as, 

a=( kL 2 -6H )/2 1 4 ,b =( 4 H-k )/L 3 , c=k/2 , 

d = e = 0 (2) 

where , 

x is the second derivative value at the end, 

H is the maximum height of the hump, 

L is half the maximum length of the hump, 

2 L is the full span of the hump. 

The equation of the other half of the hump can be derived on the 
similar lines and is found to be, 

Yz(x) = Yi(2L - x) (3) 

It must be noted that the value of « be chosen such that the 
equation should not show peak heights larger than the maximum height 
of the hump chosen. The final shape of the hump can be expressed as 
follows; 


I 


V 


Yi(X) 

Yi(2L - X) 
0 


0 S X < L 
L < X < 2L 
Otherwise 


Y( X ) 


(4) 



From the above equation it is to be understood that we are forcing 


the function to be zero outside the region [0, 2L) . Few shapes 
various x value, L and H are shown in Fig (1). 
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Fig. 1 Various hump sample profiles. Hump parameters are maximum 


for 


span, maximum height, curvature at peak, smoothness at base, 


etc. 



CHAPTER 2 


MODEL DESCRIPTION 

Next comes the question of choosing the valid vehicle model 
for simulating the results. In this chapter the various vehicle 
models along with present model, the input to the system, the 
solution methodology are discussed at length. The procedure is 
outlined for the simplest case of 1-DOF system. 


2.1 VEHICLE MODELS 


The simplest are quarter-car models which allow for only one 
dimensional, vertical (heave) motion, consisting of just the vehicle 
or sprung mass, is shown in Fig. 2. It is, moreover, a common 
practice to assume that the vehicle is moving with constant velocity 

V . 



TV 


M s 

~T 


Sprung mass. 


X .> lu c Suspension 


X 


w Jioad irregularity. 




Fig. 2 One-degree-of -freedom vehicle model . This model can b 
analyzed theoretically and was used to validate our numerical model. 
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Adding a wheel/tire sub assembly (unsprung mass) to the above 
1-DOF structure, leads to 2-DOF , quarter-car model (Fig. 3). The 
unsprung mass mode is often referred to as "wheel-hop," and is 
characterized by relatively light dampening and natural frequency 
between 8 and 12 Hz. The principal body mode, or sprung mass mode, is 
typically around 1 Hz. 


V 



1U 


M s 


K f > fcji C 1 


:n 
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t 

( ) 


Vehicle body -mass 

Suspension 
W'fuzol (a/rlo) mass 


'•-» ^ ll | TlfTK 


vy road irTKgvIa-riTiez 


(L III SI 


Fig. 3 Two-degree-of -freedom vehicle model . 


Typical half car, 2-DOF models including front and rear 
wheels, model the pitch and heave of vertical modes (Fig. 4), often 


augmented by front and rear unsprung-mass vertical dynamics. Full 





€ 



Fig. 4 Half-car 2-DOF model . This includes heave ( Z-displacement 
of vehicle centre of mass) and pitch (longitudinal rotation of vehicle 
axis ) . 


V' 



Fig. 5 Three-degree-of -f reedorn vehicle model Cincluding seat! . 
Passenger is assumed to be rigidly attached. 
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car, 3D models, on the other hand, include vehicle roll, as well as 
the pitch and heave modes (Fig. 7(a)). These models consider only 
rigid body motions, and are based on lumped-parameters, linearized 
dynamics. The suspension of a 2-DOF vehicle travelling on a randomly 
corrugated road was optimized with respect to both road holding and 
ride comfort by Dahlberg [3] . In this model the vehicle is subjected 
to stationary zero mean Gaussian random excitation. In the process of 
optimizing the elasto-damping elements of passenger car, Demic [5] 
adapted a 7-DOF vehicle model and reported a good match with the 
theoretical and experimental results. A 8-DOF model representing 
driver’s seat also was studied by Pintado [4, 13] in both frequency 
and time domains for optimizing the suspension parameters. A 2-DOF 
vehicle model with active suspension system has been used for 
optimizing suspension parameters under different conditions [4, 9, 
10, 11 and 16] . 

Several types of suspensions are used in the analysis and in 
the present work the passive suspensions is assumed. Suspension 
optimizing models are mostly linear with relatively small sprung mass 
velocities and displacement angles. 

2.2 PRESENT MODEL 

The present model assumes the following: The automobile 


model consists of a sprung mass, two independent front suspensions 
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between sprung and unsprung masses is assumed to consist of parallel 
combinations of springs and dampers. In addition, all mass assemblies 
are considered to be rigid bodies. The centers of mass of wheel 
suspension assembly is assumed at the wheel centers. In this model 
three degrees of freedom are associated with the sprung mass; bounce, 
roll and pitch. The bounce motion of all wheels relative to the 
sprung mass constitute four additional degrees of freedom. The 
vertical motion of driver’s seat relative to the sprung mass 
constitute the last and the eighth degree of freedom of the vehicle. 
Aerodynamic forces acting on the vehicle are neglected. 

The complexity of the 8-DOF model can be reduced to 5-DOF 
and 3-DOF models, for testing and validation, by shrinking the 
original 8-DOF model. For example, if we reduce the width of the 
8-DOF model to zero then it represents 5-DOF model and if length is 
also made zero, it represents a 3-DOF model. This reduction to 3-DOF 
model and 5-DOF model was shown in Fig. 5 and Fig. 6 respectively. 
The 5-DOF model is the one adapted by Essam Kassem [6] for evaluating 
dynamic effect of speed control humps. 

2.3 MATHEMATICAL FORMULATION OF 8"D0F VEHICLE MODEL 

The mathematical formulation for 3D vehicle oscillations 
model presented below (Fig. 7(a)) makes it possible to analyze the 
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following: 


a) Sprung mass vertical oscillations X3, roll w, and pitch <p 

b) Wheel vertical oscillations xi, x2, X3 and x4 

c) seat vertical oscillations xa 


The oscillatory movements of the vehicle are caused by the 
excitation generated by the irregularities of the road which the tire 
fells at the road-tire contact point denoted by wi , W 2 , wa and W4. 

Characteristics of elasto-damping elements of the vehicle 
are described by means of the parameters i = 1..4, 
j - 1..4. Referring to Fig. 7(a), the following designations apply: 


mi , mz 

ma, m4 

ms 

mo 

Ir, Ip 


li, l z, S 


-unsprung masses of the front wheel , 

-unsprung masses of the rear wheel , 

-sprung mass, 

-mass of driver’s seat including driver, 

-mass moment of inertias of the sprung mass about 
roll and pitch axes, 

-coordinates of the center of mass of the seat 
with reference to sprung mass center . 


Taking into account possible degrees of freedom in the vehicle 
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model adapted. 


the following relations are developed to define 


deformations of the elastic elements: 


Do 

— 

X5 

- 

ly. 4> 

+ 

l X 

■ w 


Di 

- 

xs 

+ 

ll.rp 

- 

s. 

w - 

XI 

D2 

= 

xs 

+ 

li. <p 

+ 

s. 

V - 

X 2 

Da 

= 

X 5 

- 

12. <p 

- 

s. 

V - 

X 3 

D4 

= 

X 5 

- 


+ 

c 

y - 

X 4 

Ds 

= 

XI 

- 

Wi 





Dc s 

= 

X2 

- 

W2 





D7 

= 

X 3 

- 

W 3 





Da 

= 

X 4 

— 

W4 






Relative velocity of the shock absorber is defined by means 
of the following expressions: 


Doo 

= 

X5 

- 

ly.'4> 

+ 

lx 

■ w 


Dp 

= 

X3 

+ 

h.<f> 

- 

S. 

V - 

XI 

Dio 

= 

X5 

+ 

U.'(p 

+ 

S. 

V - 

X 2 

Du 

= 

X5 

- 

l 2 . ’4> 

- 

s. 

V - 

X3 

Dl 2 

s 

X5 

- 

l 2. <f> 

+ 

s. 

V - 

X4 

Dia 

= 

XI 

- 

Wl 





Dl 4 

= 

X2 

- 

W2 





Dl 5 

= 

X 3 

- 

W3 





Dl<5 

- 

X 4 

— 

W4 






The following relations are developed to obtain the the 


forces in the linear elasto-damping elements'- 
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a) Tires: 

F 1 = k.21. Ds + C 21. Dl3 
Fz = kzz. D<5 + C 22. Dl4 

Fa = hza.D7 + cz3.Dis (7) 

F 4 = k.Z4. Db + C24.Di<5 

b) Suspension springs and shock absorbers: 

Fs = ten . Di + cii. Dp 
F< 5 = kl2. D2 + C 12. DiO 

F 7 = k±a. Da + ci3. Du (8) 

Fs = ki4. D4 + C 14. DiZ 

c) Driver seat’s elastic elements: 

Fo = koo . Do + Coo. Doo ( 9 ) 

On the basis of the force scheme depicted in Fig. 7(b), and 
invoking Newton’s laws, one can write the following differential 
equations describing small oscillatory motions of the vehicle model 
around the equilibrium position: 
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x o = (F o)/fflo 

XI = (F 5 - Fi)/mi. 

X 2 = (Fd - F2)//7U 

x 3 = (F? - F3)/m2 f 10 ) 

x4 = (Fa - F4)/mz 

x 5 = — (Fs + Fd + F? + Fa - Fo^/ms 

xd = (F 5 + F 7 - Fd - F 8 - Fo. Lx) S/lr 

x 7 ■ (CF? + Fa). Iz - (Fs + Fd). ii + Fo. ly)/lp 

Fi through F4 are the contact forces. When these go to zero, the car 

loses contact, where these equations are no longer valid. However 

numerically, such loss of contact is difficult to determine exactly. 
These equilibrium equations can be put in the most general form of 
the second order linear differential equations as [8, 12 and 13] , 

[M] { X( t ) } + [C]{X(t)> + LKJ{X(t)} = { F( t ) } (11) 

Where [M] , [C] and [K] are the inertia, damping and stiffness 

matrices respectively; matrices [C] and [K] might show non-linear 

• « * 

character of the dampers and springs; (X), (X) and (X) are the 

displacement, velocity and acceleration vectors in generalized 
coordinates of the multi body assemblage and the load vector (F(t)> 
is due to road irregularities, to braking forces and to centrifugal 


forces generated when the vehicle runs on a curved track. 
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The detailed description of these equations, elements of 
matrices and force vector for the case of the present vehicle model 
are presented in Appendix A. The model parameters are presented in 
Appendix B, for two vehicle types. 

2.5 SOLUTION METHODOLOGY 

2 . 5.1 TIME DOMAIN SOLUTION 

The above equations can be solved in both time and frequency 
domains using numerical methods. For getting time-domain solution 
Runge-Kutta method is used. For frequency domain solution Transfer 
Function approach is adapted. The state space representation of above 
equations yields. 


f{Xi)1 



k = 

U X2 } J 




tl] 

[K3 - LM3 A 


tC] 


{X2} 


> + 


' { 0 } 1 


C 12) 


or in the most general form. 


{X(t)}= [A] { X( t ) } + [B]{U(t)> 


C 13 ) 
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Where , 

{X} is the vector of state variables, 

{Xl} contains the generalized coordinates, both rotational and 
translational , 

{X 2 } is the first time derivative of vector {Xi}, 

* 

{Xi} is the first time derivative of vector {Xi}, 

* 

{Xz} is the second time derivative of vector {Xi}, 

{U(t)> is the input vector. 


2 . 5.2 FREQUENCY DOMAIN SOLUTION 


The frequency domain solution can be obtained by following 
state space approach as follows. The solution to (13) can be assumed 
as 


(Y(t)> » [A] { X( t ) } + [D] { U( t ) } 


( 14) 


Taking Laplace Transforms of both equations (13, 14 ) we get. 


s { X( s ) } - { X( 0 ) } = [A] { X( s ) } + [B] { U( s ) } 
(Y(s)} = [C] { X( s ) } 


with zero initial conditions i.e., { X( 0 ) } = {0}, we get 
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s (X(s)> = [A] { X( s ) } + [B] ( U( s ) } 

(Y(s)> = [CJ { X( s ) } 

solving for (X(s)}, the above equation give 

(X(s)} = (s [I] - LAjVlBJ (U(s)l 

and 

(Y(s)> = [C] ( $ [I] - [A]) 1 [BJ { U( s ) } 


or , 

{Y(s)> = [ [G( s )] {U( s )}] (15) 

where , 

[G(s)]= [ [C] ( s [I] - [A] )* [B] ] (16) 


is the Transfer Function matrix defined as the ratio of Laplace 
Transform of output to that of input. And, 

s is a complex quntity and is equal to jo> 

where j - i -1 , and co is the angular frequency expressed 
in rad/sec, <0 = 2 nf , / is the frequency in Hz, 

[I] is the unit matrix, 

[0] is the zero matrix, 

(U(s)> is the Laplace Transform of input i • » no (U(t)) 

(Y(s)} is the Laplace Transform of output i.e., £(Y(t)) 

{X(s)> is the Laplace Transform of X( t ) 

[Cl is a matrix, in the present case it is unity. 
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Because l C3 is unity, (X(s)} = (Y(s)}. 

This transfer function, as seen, is a constant function of s 
for the given system and therefore it is a constant function of a as 
long as the system is time invariant i.e., the system parameters does 
not change with time. It is independent of input to the system. The 
output of the system coordinates in frequency domain can be obtained 
from (15). 

2 . 5.3 TIME DOMAIN INPUT 

Assuming the vehicle velocity is not affected by the road 
geometry, front and rear tire inputs (( wi, W2 ) and (w3, W4 )) are 
related by the following expressions: 


wi, wz = Yl(t) (17) 

and , 

W3, w4 = Yz( t — t ) (18) 


Where the time delay r depends on the the wheel base and vehicle 
running velocity (K). 

Wheel-base W 


I 


V 


V 


( 19) 
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Therefore the input to the front and rear tires in time domain can 
be expressed as. 


wi = Y(t) 


Yu = A(t-to) 4 + B(t-to) 3 + C(t-to)^ 
- Y 12 = v4(ti-t) 4 + S(ti-t) 3 + C(ti-t) 2 , 
Yl3 = 0 


to < t < L/V + to 
L/V < t < 2 L/V 

Otherwise 


and 


wa = Y(t - t) 


f Y 21 
Y 22 
Y 23 


A(t-t2) 4 + B( t“t2 ) 3 

/4(t3-t) 4 + S(t3-t) 3 

O 


+ C( t~t2 ) Z , 
+ C(t3“t)f 


W < t < ( WtL )/V 
( W+L )/V < t < ( W+2L )/V 
Otherwise 


( 20 ) 


Where , . 

4 3 _ Z 

A =■ a V , B = b V , C = c V . (a, b, c are defined in (2)), 
W = It + I 2 is the wheel-base of the vehicle, 
to > 0 , is the time taken by the car to reach the bump, 
tl = to + 2 L/V, t2 — to + W/K, t3 — to + ( W + 2 L)/V, 

V is the forward velocity of the vehicle. 


2 . 5.4 FREQUENCY DOMAIN INPUT 



To calculate the input vector (U(t)> in frequency or s 
domain, we need to consider the Laplace transform of { U( t ) } and here 
for the inputs wi, W2, W3 , W4 the Laplace transforms are as follows: 
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and , 


ex 

WKS), W 2 iS) = £(Y(t)) = J Y(t).e -St dt 


( to+L/V > 

Yi*( t ) .e st "dt 

to 


+ 


o 

(2L/ V) 

Yi2( t ) .& -L dt 

L/V 


cx 

W3<S>, W4<S> = £,( Y( t-r ) ) = JY( t-r ) .e ~ S kt 

o 


< t2+L/V) 

Yii(t).e" sL dt 

12 


+ 


(l 2+2L/V) 

Yi2( t ) .e Sl dt 

< t2 + L /V > 


(21 ) 


( 22 ) 


from first shifting property of Laplace transforms, W3 can be 
expressed as, 

W3<s>, W4<s) = e ST .£( Y(t)) = e ° T „wi<s> (23) 

2.6 TRANSFER FUNCTION FOR 1“D0F VEHICLE MODEL CASE 


The transfer function approach for obtaining frequency 
domain solution for the case of 1-DOF vehicle is explained below. The 
governing dynamic equilibrium equation for the vehicle model shown in 
Fig. 2 can be written as. 


m. z* + c z + tez 


k V) + c w 


( 24) 
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Where , 

m - Ms is the sprung mass (Kg), 

k = Ki is the suspension spring constant ( N/m ) , 

c = Ci is the suspension damping rate ( N-sec /m ) , 

z, £, z are displacement, velocity and acceleration of mass m, 

( » dot indicates derivative with respect to time ), 

li' indicates the input excitation to the system (road 
irregularities ) . 


Transforming this equation into state space representation, the 
matrices in (12) represents. 


' o 

• 

1 

[B] = 

o 

X 

0 

and (U(t)> = 

'1 

-k/m 

-c/m 

J 

2X2 

h/m 

c/m 

2x2 



T 

and the state vector {X} is [ z z'] 

Zxi 


T.F = [G] = Is. I - 


'k/m 

c/m 

1 

Gil 

Gi2 

s. k/m 

s, c/m 

A 

G2i 

GZ2 


{Y(sO) = [G] (U(«)> = 


' h/m. M.s) +c/m v(&) 
sCk/m s )+ c/m itXs)) 


1 _ 

A 
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where A is the determinant of [G] and is equal to. 


A 


2 


S 


+ s, c/m + k/m. 


Displacement of mass m. 


z(s)= Yi( s ) = (k/m..nKs) + c/m. w( s ) ) / A 


If w( t )= Ao sin(w.t) then, 


n<s) = £(ii<t))= ft( sin( a>.t ) ) = — 1 — 

( « + s ) 

and 

= £(w(t)) = £(w.cos(oo.t )) = — A ^' a) '~ -- 

(to + s Z ) 


then z(s)is given by. 


z( s) 


Wm. (M + c. sJVCw 2 + 


A 


The transfer functions Gu, G 12 for this system case are shown in 
Fig. 8(a) and the corresponding input in s domain is shown in Fig. 
8(b) and the resulting output in s domain is given in Fig. 8(c). 
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CHAPTER III 

HUMP SELECTION - A DISCUSSION 

Now the question is "What are the suitable parameters of a 
hump that best meets the objectives, if the limiting speed and the 
vehicle type are specified?" Considering the conflicting factors viz. 
road-holding and ride comfort, as evaluated by conducting tests on 
level roads, it is not possible to directly use them for case of 
travel across a hump in evaluating the acceptability of the hump. As 
no test data is available regarding these factors, in this work, the 
dependency of the dynamic response of vehicle and driver on road 
sections with hump is presented. 

In this chapter the general criteria for evaluating the 
performance of vehicle and driver are discussed. The criterions 
modeled in this thesis are described. Simulation results are presented 
in the form of graphs for various travel considerations of vehicle 
across different humps. 

3.1 ROAD HOLDING PROPERTY 

Optimal road holding can be seen as the minimal probability 
that the magnitude of road-wheel contact force will be larger than a 
given level. When a vehicle crosses a hump at a speed, for safety, 
the tire should not lose contact with the ground. This phenomena of 
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contact loosing can be quantified best with the road - wheel contact 
force (7). It can be seen that these forces (7) are proportional to 
the tire vertical displacement and velocity. Under severe cases, the 


wheels can 

be 

deflected 

by large amounts 

causing 

damage 

suspension . 

The 

contact 

forces F, ( L=1 ,4 ) in 

i 

equation 

(7) 


deviations from the nominal contact force due to the weight of the 
vehicle. Thus, even if F, is negative, so long as (F. + P ) is 

L L L 

positive, contact is not lost. P, is the total load share on any tire 

i 

(For example in Fig. 19, assuming P,= W/4 on each tire, the contact 

(. <■ 

force is negative only marginally). 

To select a hump from this criteria and result given in this 
work, one have to have a limiting load on the tire. If at all one 
wishes to select a hump from these results, he has to decide himself 
the allowable/acceptable limits of these factors. For example, if we 
limit ourself that the wheel-road contact force should not be more 
than 1.75 times that in static condition when the vehicle crosses the 
hump, for the vehicle Make-1 presented in Appendix B, then the 
assumed velocity at which contact is lost can be seen to be 6 m/sec 
(21.5-KM PH) (Figures 37 and 38). 


i 

i 
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3.2 HUMAN COMFORT 

With human comfort point of view, an ideal hump should be 
such that at and below the limiting speeds ride is to be comfortable 
and vice-versa. There are many a different ways to quantify the ride 
quality. The simplest one is based on the calculation of the 
root-mean-square( rms ) value of the vertical acceleration typically 
measured at the driver’s or passenger’s seat location. A useful early 
study in this regard was published in 1978 [14] . The authors draws 
conclusions based on field study, states that "the values of these 
magnitude weighted rms values will range roughly from 0 to 0.04# for 
smoothC interstate highway) rides, 0.04 to 0.06# for medium rides, and 
above 0.06# for rough rides." 

The rms quantifiers do not reflect the frequency dependence 
of human sensitivity to vibrations. To account for this aspect of 
ride comfort, the International Standard Organization (ISO) has 
developed a standard , which considers the duration of human exposure 
to the vibrations (ISO 2631). According to the above standard, the 
region of greatest human sensitivity to vertical vibrations lies 
between 4 and 8 Hz, which roughly include various resonance 
frequencies of human internal organs. The above measures are 
primarily based on decoupling the degrees of freedom typically 
performed in only one direction at a time. It is not clear the 
sustained input models such as these would apply to non periodic and 
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shock loading conditions as in a road hump. To select a hump from 
this criteria and simulation results presented here, one have to have 
a limiting/tolerable acceleration limits which can be directly 
comparable. Here, frequency dependency of vertical acceleration and 
displacement at driver’s seat location are also simulated because, if 
at all any text is going to be conducted for this case, it would be 
easy to measure and compare these quantities. 

3.3 OTHER ASPECTS 

One often introduces a constraint in the form of rms 
suspension stroke or so called rattle space. Rattle space is a term 
used to indicate the clearance allowed for the vertical movement of 
wheel. High rattle space indicates the possible touching of the axle 
with the chassis or the car body, causing possible damage. This 
aspect is important for case of travel across a hump because, there 
may be a possible large and sudden deflection of suspension which may 
damage vehicle itself. This aspect is very significant for designing 
suspension, here it is only considered to study the variation of the 
packaging constraints. 

3.4 PRESENT MODEL ACCEPTANCE CRITERIA 

In the present work the following performance criteria has 
been modeled which are being widely adapted for various purposes of 
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evaluating the vehicle performance [1, 2, 3, 4, 6, 9, 10, 11, 13, 14, 
15 and 16]. In time domain, 

a) . Vertical acceleration at drivers’s seat (related to ride 
comfort ) . 

b) . Suspension working stroke or rattle space (related to 
packaging constraints of the vehicle body). 

c) . Road-wheel contact force variation (related to road holding 

property ) . 

and in frequency domain, 

d) . Vertical acceleration of the driver . 

e) . Vertical displacement of the driver. 

From Fig. 9, it is clear that as vehicle velocity is 
exceeding certain value, the rms vertical acceleration of driver is 
decreasing indicating a good drive comfort. This encourages the 
motorists to travel at much higher speeds and the hump in question is 
becoming a overhead. This figure also compares two types of vehicles, 
one with larger mass (Make-2) and the other with lesser mass 
(Make-1). It is clear that the vehicle with larger mass is having 
much higher dynamic effects at a given speed compared to the lesser 


weight one. 



uol ■^•' E J °1 
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Fig. 11 Effect of crossing speed on the rms rattle at front end. 
CMake-1 and Make-2 refer to vehicle data). 



Fig. 12 Effect of crossing speed on the rms rattle at rear end. 
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Fig. 13 Effect 

of crossing speed on the rms 

road 

wheel contact 

force 

at front end. 
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Fig i.4 - Effect of crossing speed on the rms road wheel contact force 
at rear end. , C Make-1 and Make-2 refer to vehicle data}. 
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Fig. 10 displays rms vertical displacement of the driver 
This displacement is also decreasing beyond certain velocity level 
Figures 11 and 12 show that the rms rattle space is also decreasinj 
with increasing speed, and this magnitude is higher for lesser mass 
car than the other. On the other hand, the rms wheel-road contacl 
forces show an increasing effect with velocity. That means this 
contact force may cross the permissible limit and thus may nol 
satisfy road-holding criterion. Often it may be possible from thes< 
plots to draw the fly-off velocity assuming assuming the limit valu< 
of the contact force. Moreover, large mass vehicles experience 
relatively large contact forces and is very likely to loose contacl 
at lower velocities compared to lesser mass vehicles. 

3.5 EFFECT OF VARYING THE WIDTH OF THE HUMP 

Five different hump widths are examined. These are: l.o 
1.25, 1.5, 1.75 and 2.0 m. Here hump height(0.15 m) and crossin 
velocity(5.0 m/s) are kept constant. Figures 15-20 show the effec 
of varying the width of the hump on the various aspects discussed i 
the previous section. A general conclusion can be drawn upo 
examining the response figures which say that the wider the hump th 
more is the dynamic effect on driver and the rattle space. Where a 
tensile part of the road holding force is decreasing making th 


drive safe. 
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Fig. 19 Effect of varying the width of the hump on the wheel-road 
contact force at front end. Road contact force approaches to a 
constant value as hump width is increasing. 3hCSC contact 4'ofCCS Ofe 
deviations from the. nominal contact -force, due to tht weigte of the. vehicle 
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Fig. 20 Effect of varying the width of the hump on the wheel-road 
contact force at rear end. 
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Fig. 21 Frequency response of driver’s vertical acceleration for 
various hump widths. In general, at higher frequencies the same 
effects are less bearable. 
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comfort levels of human are non linearly reducing functions of 
frequencies. These do not apply to shock loads, but the general 
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The frequency domain magnitudes of vertical acceleration and 
displacement of driver (Fig. 21 and 22) also show an increasing 
effect with increase of hump width. This is because of the increasing 
crossing time with increasing width (i.e., the disturbance from road 
is acting for a longer time). 

3.6 EFFECT OF VARYING HEIGHT OF THE HUMP 

Five different maximum heights of humps are examined. These 
are 0.1, 0.15, 0.2, 0.25 and 0.3 m. Figures 23-29 show the effect of 
varying the hump height in time domain and figures 30 and 31 show 
this effect in frequency domain. As in the previous case the effect 
is very much the same. 

3.7 EFFECT OF CROSSING SPEED 

The significance of high speed is that the vehicle crosses 
the hump in shorter time. This can be noticed in all figures (32 
38). Five different constant crossing speed are tried. They are 2, 3, 

4, 5 and 6 m/sec. All aspects considered so far are having a tendency 
to decrease with increasing velocity except the road-tire contact 
force . 

The tendency of increasing road-tire contact force is 
helpful in deciding whether a particular speed is allowable or not 
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Fig. 27 Effect of varying the height of the hump on the 
wheel -road contact force at front end. These contact forces are devi- 
attans ■from nominal contact forces due tt the. weight of the. vehicle- 



Fig. 28 Effect of varying the height of the hump on the 
wheel -road contact force at rear end. 
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Fig. 29 Frequency response of driver’s vertical acceleration for 
various hump heights. 
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This result indicates that at speeds causing forcing frequencies 
greater than the natural frequency of the system, less energy is 
transferred to the system. Since permissible discomfort usually 
Increases with acceleration, these results indicate greater comfort at 
higher speeds . 
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Fig. 39 Effect of varying the crossing speed Dn the wheel . road 
contact force at front end. These contact Forces are dLea.'/iationS from 
the nominal contact force <io.c to tna weight of the vehicle. 
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Fig. 30 Effect of varying the crossing speed on the wheel-road 
contact force at rear end. 




Fig. 37 Frequency response of driver's vertical acceleration for 
various crossing speeds. 



Fig. 38 Frequency response 
various crossing speeds. 


of driver's vertical displacement for 


from safety point of view. It will be helpful to explain the 
decreasing effect of acceleration in frequency domain (Fig. 39). From 
this figure it is clear that for low velocities the magnitude of 
acceleration is dominant at low frequencies only. This can be 
tolerable because of low frequency. On the other hand, this magnitude 
is dominant at high frequency for high velocities, which is 


intolerable . 



CHAPTER 4 


CONCLUSIONS AND SCOPE FOR FURTHER WORK 


We started the hump design process with the objective of: 
having little or no effect at low speeds and significant effect at 
high speeds. The following conclusions can be drawn based on the 
simulation results. 

(1) The dynamic effect of humps vary considerably with the type 
of vehicle. The vehicles with large mass undergo comparatively 
greater dynamic effects. 

(2) The hump parameters (height and width) also affect the 
dynamic response. 

(3) For a given hump height (hump width) and crossing speed 
discomfort and suspension rattle space increases as hump width (hump 


height) increases. 
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(4) For a given hump height (hump width) and crossing speed, road 
contact force approaches a constant value as hump width (hump height) 
increases. This constant value is the same as that due to the 
stationary vehicle weight. 

(5) Vertical acceleration of driver decreases as speed increases. 
Ride comfort is dependent on the vertical acceleration of the driver; 
the higher the acceleration, the less the ride comfort. Hence, ride 
comfort is likely to be better as speed increases. 

The rms value of vertical acceleration increases with speed up 
to a maximum value (which depends on the vehicle type); as speed 
increases further the rms acceleration decreases. Intuitively this 
appears anomalous, but less energy is transferred to the vehicle as 
the frequencies increase well beyond the natural frequencies of the 
system. This is the inverse of the behavior a speed control device 

should have . 

(6) Magnitude of road-holding force is directly proportional to 
the velocity of the vehicle. This is larger for vehicles with larger 


mass (Fig. 13). 
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(7) The magnitude of road holding force increases with speed up 
to a maximum value (which depend on the vehicle type) indicating good 
road holding, as speed increases further the magnitude of road 
holding force increases largely (thus increasing the hazards) while 
crossing a hump. 

(8) Conclusions 4 and 6 further reveal that in general the 
comfort increases and. safety decreases with increasing crossing speed 
of the vehicle. 

(9) Given the limiting/allowable contact forces these simulation 
results give the fly-off velocity of vehicle for travel a cross ahump. 

(10) Frequency dependency of driver’s vertical acceleration and 
displacement are also presented for a class of humps. These can be 
directly used to test the acceptability of a hump if any test data in 
this field is available. 

What do these results imply for hump design. Specifically 
it appears that rider comfort cannot be a factor in hump design, 
since comfort improves at high speeds. On the other hand, rattle 
space and road holding (i.e., safety) decreases as speed increases. 
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These results appear to be contrary to normal intuition, and 
experience while travelling over humps. A possible explanation to 
these contrary results evolved after the work was finished. The 
reasons can be explained with reference to a 1-DOF vehicle model 
(Fig- 2). The following figure shows the frequency domain input for 
two different crossing speeds of the vehicle and the transfer 
function. Typically, at higher velocity, frequency domain input has a 
wider frequency spread, but lesser amplitudes. 



Figure showing the frequency domain input for two different 
crossing speeds of the vehicle. 


Magni lude 
o f 

Tran sfer Function 


Figure 



showing transfer function for 1-DOF vehicle model. 


It can be seen from these figures, the input (hump) is such 
that at velocities Vi and Vz, both are present even after the natural 
frequency of the system, except that their magnitudes are different. 

The vehicle while travelling the at a lesser speed (W> absorbs more 

> 

1 
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energy compared to travel at a higher velocity velocity ( Vz ) since 
the amplitude of frequency domain input is higher. This could be the 
reason why the results are contrary to intuition. 

On the other hand if we choose a hump profile whose 
frequency domain input for two different velocities appears like the 
one shown in the below figure, then it would be possible to achieve 
the objective of design of a hump. 



Figure showing the frequency domain input for two different 
crossing speeds of the vehicle. 

Further work needs to try to choose humps such that at low 
speeds the forcing input appears as the input A and at higher 
velocity as the input B. In this type of hump profile the true 
objective of hump design may be achieved. 

Another correction that would be desirable in future work 
could be to prevent the total contact force from becoming negative. 
This would imply that the tire can be pulled down by the ground which 
is impossible. To correct for this, the forcing function should be 
nonlinear in the sense that the function cannot go below zero. 
However this would lead to considerable complexity of since one would 
then need to model the effect of impact after loosing contact. 


APPENDIX A 
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STATE-SPACE EQUATION 


From the car body, front and rear suspensions, tires and 
seat, the state-space equation of the vehicle modal in matrix form is: 

X(t) = A.X(t) + B.U(t) 

The 16-components of state variable vector is: 

X = [ Xo, XI, X2, X 3, X4, X5, >//, (p, Xo, XI, X2, X3, X4, X3, $>] T 

The input vector is: 


U = [wi, Wl, W 2 , W 2 , wa, W3, W4, W4, 0, cf t p, 0] 

where cf is the force due to centrifugal effect and is equal to 
2 

(mis.K R/h ) and p is the force that arises due to acceleration and 
deceleration effects and is equal to (ms.a.h). Here R is the 
reciprocal of curvature of turning circle of the vehicle, h is the 
height of the center of mass of the vehicle sprung mass, and a is the 
acceleration or deceleration m/sec ) of the vehicle. Here uniform 
velocity is assumed (c / = p = 0). 


The state equation matrices can be expressed from (12) and (13) as. 


A * 



- 


■ * 

0 

I 


Bn 

8X8 

8x8 

B = 


[ -M" 1 . K] 

*• fl X 0 

[ -M _i . C] 

1 J 0X0 


B 22 




“ * 


1(5x8 



where 
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B “ ' 101 SXB 822 


bl b2 0 o o 

0 o ba b4 o 

0 0 0 o bs 

ooooo 

ooooo 

ooooo 

ooooo 

ooooo 


0 0 0 

ooo 

bd 0 O 
0 b? be 

ooo 

ooo 

ooo 


0 


0 


0 


8 x0 


Ma br .i x M - It is a diagonal mass matrix with elements: 


mii = mi, m 22 = mz, m33 = ma, m44 = nu, 
mss = ms, m<3<3 = Ir, m?? = Ip, mae = mo. 

Matrix K: It is a symmetric matrix with spring constants and its 
elements are as follows: 


ki,i a -(Mu + Mzi ), ki,s = Mu, ki ,<5 = -S. Mu, ki ,7 = Mu. li 
k 2,2 = -( M12 + M 22 ), k 2 ,s = Miz, k 2,<5 = S.kiZ, k 2,7 = Miz. li 

ka,a = -(Mia + M 23 ), ka,s = M 13 , ka,d = -S. Mia, k 3,7 = -M13. Iz 

k4,4 = -( Ml4 + M 24 ), k 4,5 = Ml 4 , k4,(5 = S.k 14, k4,7 = -Mi4. 12 

ks,s = -( Mu + Mi 2 + Mi 3 + Mi 4 + Mo ) 

ks ,<5 a S. (Mu - M 12 + Ml3 - Ml 4 ) + lx. ko 
ks,7 = -il. ( Mu + Ml 2 )+ I2. ( Ml 3 + M.14 ) - Mo. ly 
ks,a = Mo 

<y 2 

L.^ ^ ~ _<£**/ _i. Uj, 4. a + kiA. ) - lx. Mo 


57 

k<5,7 = -S. lt,( ki 2 - kii ) - s. L 2 . 0.13 - ki 4 ) + lx. Iv.kc, 
k<5,Q = - Ik. ko 

l<7,7 = “ll. ( &11+&12) - Iz. (kia+ki4) - ly.ke 
k7,B = ly. ko 
ka,8 = -ko 

Matrix C- It is a symmetric matrix with damper coefficients and its 
elements are as follows: 

C 1,1 = -(Cii + Czi ), Ci, 5 = Cii, Ci, a = -S. Cn, Ci , 7 = Cu. li 
CZ,2 = “(Cl2 + C22), C2,5 = Cl2, C2,<3 = S.C 12, C2,7 = ClZ. li 

C3,3 = “(Cl3 + C23), C3,5 = Ci3, C3,<5 = “S. Ci3 , C3,7 = “Cl3. Iz 

C4,4 = “( Cl4 + Cz4 ), C4,5 = Ci4, C4,<5 = S. Ci4 , C4,7 = -Ci 4 . Iz 

C5,5 = -( Cll + Cl2 + Cl3 + Cl4 + Co) 

C5,<5 = S. (Cii ~ CiZ + Cl3 “ Cl4 )+ ly.. Co 
CS,7 - + Ci2)+ 1 2 . ( Cl3 + Cl4 )- Co. Ly 

C5,a = Co 

2 2 

C<S,<5 = ~S. ( Cll + C 12 + Cl3 + Cl4 )~ lx. Co 

C<3,7 = “S. Ll. (ClZ - Cll )~ S. 12. (Cl3 - Cl4 )+ lx. ly. Co 

C<5,0 = — lx, Co 

C7,7 = “ll. (Cn+Cl2)- ll, (C13+C14 )“ ly.Co 
C?,B = ly. Co 
Ca,a =: “Co 

WOT E s I_r> the abo.e matrices element specified other than symmetric elements are all zeros. 
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Data for two vehicle models (Make-1 and Make-2) used in 


present work are given in the below Table 


Parame t er 

Symbol 

Uni ts 

Make-1 

Make-2 

Car body mass 

m s 

Kg 

1 090 

2 000 

Car body mass moment 





of Inertias: 





In Pitching 

Ip 

2 

Kg~m 

1 680 

3 200 

In Rolling 

Ir 

2 

Kg-m 

250 

510 

Front unsprung mass 

mi ,m 2 

Kg 

27.5 

25.0 

Rear unsprung mass 

ma , m4 

Kg 

25.0 

25.0 

Suspension stiffness: 


N/m 



Front 

ki l ,kiz 


7 500 

20 000 

Rear 

kz i , kzz 


7 500 

20 000 

Suspension damping: 


N .s/m 



Front 

C 11 ,C 12 


1 000 

1 500 

Rear 

CZ 1 ,C 22 


1 000 

1 500 

Tire stiffness 


N/m 



Front 

kz i , kz 2 


75 000 

90 000 

Rear 

kz 3 , kz 4 


75 000 

90 000 

Tire damping 


N .s/m 



Front 

C 2 1 ,C22 


200 

275 

Rear 

C23 ,C 24 


200 

275 

Wheel-base 

w 

m 

2.5 

2.5 

Wheel-track 

s 

m 

0.7 

0.7 

Front length 

li 

m 

1 .4 

1 .3 

Rear length 

Iz 

m 

1 .1 

1.2 

Seat total mass 

! mo 

Kg 

100 

100 

Seat spring stiffness 

ko 

N/m 

1 500 

10 000 

Seat damping coeff. 

Co 

N .s/m 

150 

500 

Seat X-position 

lx 

m 

0.4 

0.5 

Seat Y-position 

ly 

m 

0.35 

0.35 


the 
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